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Abstract. We demonstrate that the thermodynamics of one-dimensional Lieb- 
Linigcr bosons can be accurately calculated in analytic fashion using the polylog 
function in the framework of the thermodynamic Bethe ansatz. The approach does 
away with the need to numerically solve the thermodynamic Bethe ansatz (Yang- 
Yang) equation. The expression for the equation of state allows the exploration of 
Tomonaga-Luttinger liquid physics and quantum criticality in an archetypical quantum 
system. In particular, the low-temperature phase diagram is obtained, along with 
the scaling functions for the density and compressibility. It has been shown recently 
by Guan and Ho (arXiv:1010.1301) that such scaling can be used to map out the 
criticality of ultracold fermionic atoms in experiments. We show here how to map 
out quantum criticality for Lieb-Liniger bosons. More generally the polylog function 
formalism can be applied to a wide range of Bethe ansatz integrable quantum many- 
body systems which are currently of theoretical and experimental interest, such as 
strongly interacting multi-component fermions, spinor bosons and mixtures of bosons 
and fermions. 
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1. Introduction 

The one-dimensional Lieb-Liniger model [1] of interacting bosons is arguably the 
simplest Bethe ansatz integrable model [2j |3]. It is one of the most extensively studied 
many-body models of cold atoms|3 Experimental studies of cold bosonic atoms over a 
wide range of tunable interaction strength between atoms are in agreement with results 
obtained from the integrable model [8j [9j [10]. Moreover, the temperature dependence 
of this system provides a further experimental test [11] of integrable theory in the form 
of the Yang- Yang equation [12] in the weak coupling regime. 

| See, for example, references [H El [7]. 
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In principle, the Yang- Yang method gives exact results for the thermodynamic 
properties of integrable many-body systems at finite temperature [13]. However, 
depending on the particular model under investigation, the Yang- Yang method involves 
either a finite or an infinite number of coupled nonlinear integral equations, which 
hinders access to thermodynamic quantities from both the analytical and numerical 
points of view. It is a formidable task to extract exact low temperature results for 
many-body systems of this kind. 

Here we further develop a polylog function method to derive the equation of 
state of homogeneous degenerate quantum gases of cold atoms in the framework 
of the TBA formalism. The polylog function method has been recently applied 
to the one-dimensional attractive Fermi gases of ultracold atoms [IJJ [151 HSl fTT] . 
The thermodynamics and full phase diagrams can be calculated analytically via this 
approach. For example, for spin- 1/2 attractive fermions, the universal Tomonaga- 
Luttinger liquid behaviour is readily identified from the pressure given in terms of 
polylog functions [14] . Here we investigate the thermodynamics of the Lieb-Liniger 
model for spinless bosons. The high precision of the equation of state obtained for 
strong interaction provides exact physical properties of this model at zero and finite 
temperatures. Following [18j [H] , the universal Tomonaga-Luttinger liquid physics [19] 
follows with the help of Sommerfeld expansion in the low temperature limit. It was 
recently shown [17] that dynamical critical exponents and correlation exponents can be 
read off the universal scaling functions for thermodynamic properties of one-dimensional 
fermions. As we shall see below, the polylog function method is also accurate, simple and 
convenient in the study of thermodynamics and quantum criticality of one-dimensional 
bosons. 

2. Thermodynamics of one-dimensional bosons 

The Lieb-Liniger model [1] is described by the Hamiltonian 



in which N spinless bosons, each of mass m, are constrained by periodic boundary 
conditions on a line of length L and #id = h 2 c/m is an effective one-dimensional coupling 
constant with scattering strength c. We will be particularly interested in the low density 
case, i.e. the dimensionless parameter 7 = c/n is finitely strong, where n = N/L is the 
linear density. 

It was for this model that Yang and Yang [12] introduced the particle hole 
ensemble to describe the thermodynamics of the model in equilibrium, which is later 
called the thermodynamics Bethe ansatz (TBA), see [13]. In terms of the dressed 
energy e(k) = Tln(p h (k)/ p(k)) defined with respect to the quasimomentum k at finite 
temperature T, the Yang- Yang equation is 




(1) 
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where p is the chemical potential and e°(k) = -^k 2 is the bare dispersion. The dressed 
energy e(k) plays the role of excitation energy measured from the energy level e(kp) = 0, 
where kp is the Fermi-like momentum. The equilibrium states are described by the 
equilibrium particle and hole densities p(k) and p h (k) of the charge degrees of freedom, 
which are subject to the equation 

p( k ) + = ir + ir r ^ dq il7! q) v W 

2tt 2tt J_ oc c 2 + {k- q) 2 
The pressure p(T) and free energy F(T) are given in terms of the dressed energy 



by 

n(T) = 

2tt 



p(T) =— I dk ln(l + e~ e{k)/T ) (4) 



rp poo 

F(T) = pn / dk ln(l + e~ e{k)/T ). (5) 

2vr J_ 00 

It is important to observe that the log term in the Yang- Yang equation (j2]) vanishes 
exponentially for a large quasimomentum k > kp at low temperatures. For T — > 0, the 
ground state properties are determined by the so-called dressed energy equation 

= fr^f^. (6) 

The integration boundary Q = kp at zero temperature, i.e. e(±Q) = 0. The pressure is 
given by 

i r Q 

p=-— / dke(k). (7) 



2vr. 

This observation naturally suggests an expansion in powers of 1/c in the TBA (JSJ), i.e 



e(k) = e°(k)-p+C—^- f dqe{q) 
irc z + k 2 J Q 

By way of a straightforward calculation with iterations through the pressure (J7|) and 
the condition e(±Q) = 0, we find the relation 



which gives 



1 4 12 (^tt 2 -32) . 
E w -n 3 n 2 1 - - + — + / (10) 

5 sy z 1 



as a highly accurate expansion for the ground state energy per unit length (in units of 
h 2 /2m). 

For strong repulsion, Lieb-Liniger bosons behave like ideal particles obeying 
fractional statistics |20j. The gas reaches the Tonks-Girardeau regime with two-particle 
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local correlation g^> = (x) 2 ^/ (x) 2 ) / n 2 — > 0. Using the Hellmann-Feynman theorem, 
to leading orders the zero temperature two-particle local correlation follows as 

which agrees with the result obtained from form factors of the Sine-Gordon model in 
the non-relativistic limit [2T] . 

Turning to finite temperatures, the Yang- Yang equation can be expanded to 0(l/c 4 ) 

as 



e(k) = £ ° W - „ - fSl + _X_ £ dq f ln(1 + e - 



where 



and 



«°(*) - A* - " ^^r Li f (-^° /T ) ( 12 ) 

C + K 2y/TlC 6 (h?_\2 2 

V2m/ 



2p(T) V*/ 2 

Ao = /i + ^-i57^ (13) 



l ^) = Ef (14) 

fc=l 

is the standard polylogarithm function [22]. From the dressed energy (TP2"]) . the pressure 
at finite temperatures follows in terms of the polylog function as 



p(T) "-^2^ T§Li l ( - eA/T) 




Li3(-e A / T ) 



(15) 



where now 



A = a* + ^ + n * r " 3 Li5(-e^ r ). (16) 

\2m) 



The result ( 1151) is essentially a high precision equation of state for Lieb-Liniger bosons. 
Figure 1 and Figure 2 show the excellent agreement between the analytic result ffT5l) 
and the result obtained by numerically solving the TBA (J2]) at finite temperatures. 

We touch briefly now on some physics of the model. For the grand canonical 
ensemble, the two-particle local correlation can be obtained from the free energy per 
particle / = F{T)/n as g^ = jgr (df/dj) n T [231 HBJ. The results also allow exploration 
of the crossover from the universal Luttinger liquid regime to the decoherent regime 
where the linear dispersion in the low-lying excitations is destroyed. Considering the 
low temperature limit, i.e. <C 1, the pressure per unit length with fixed n is, 

to leading orders in T, given by 

(17) 
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which follows directly by Sommerfeld expansion [IB] • Here we have ignored higher order 
corrections than 1/c in the temperature-dependent terms. In the above result p Q and 
Hq are the pressure and chemical potential at zero temperature (in units of h 2 /2m) 

2 3 , / 6 24 (fvr 2 -80) 
Po ~ ^n 3 n 2 1 - - + — + ±1 '- 

9 o / 16 20 /64 9 \ 1\ 
w , H „ V ^__ + _ + ^_ 6 4j_j. (19) 

Moreover, to leading order, the free energy follows as 

F(T) „ Eo _ ^ni m 

bnv c 



where the central charge C = 1, Eq is the ground state energy ( ITOl) and t> c ~ 
^1 — i + i| J is the sound velocity. The result (|20|) is as expected from conformal 
field theory arguments for a critical system, i.e. for a system with massless excitations. 
This implies that for temperatures below a crossover value T*, the low-lying excitations 
have a linear relativistic dispersion relation, i.e. of the form u(k) = v c (k — kp). If 
the temperature exceeds this crossover value, the excitations involve free quasiparticles 
with nonrelativistic dispersion. This crossover temperature can be identified from 
the breakdown of linear temperature-dependent specific heat. Figure 2 indicates this 
universal crossover at a temperature T* ~ 1 (in units of ks) from a relativistic dispersion 
relation to a nonrelativistic dispersion. 

For this simplest of models, the chemical potential drives a quantum phase 
transition from a vacuum phase into a Tomonaga-Luttinger liquid phase at zero 
temperature. We now further refine the equation of state to map out the universal 
low temperature quantum phase diagram. 



3. Equation of state, scaling functions and phase diagram 



For convenience in deriving the equation of state, we introduce an energy scale Eq = 
h 2 c 2 /(2m). Equations (EES} and flE) can then be written in terms of the dimensionless 
quantities Ji = h/eq and T = T/e as 



E C 



fl 



— Li3 i — e 

2^ ^ 



1 + 



1 2 



— =Li3 ( 
2v^ 5 



and 



A = p,- 



Li: 



Jo/T 



71 



) + 



lLhu(-e A °/ f ) 



20F 



(21) 



(22) 



with A = fi — ^=Li| (— d x ' T ). After a straightforward calculation, the density and 
compressibility are given by 



n« -^Lhf-e^) 



2^ 



Ul (-e Ao / T ) + -Lh(-e A °/ T y 
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7T2 



2Jtt 5 



(23) 



2^0^ \/T 



Li_ i ( — e 

2 



3T5 t 

-Lii (— e 

'7T 2 



+ ^ Lil (_ e ^/T) 

7T 2 



(24) 



In this model there exists one critical point at fi = /i c = 0, i.e. a quantum phase 
transition from the vacuum phase into a Tomonaga-Luttinger liquid occurs at /i = \i c 
at zero temperature. Near the critical point fi c we find the density obeys the universal 
scaling form [21] 



n(T,fi) - n (T,/i) « Ti 
for which the scaling function is T{x) 



-l— 



r [ } UZ 

=Lh(-e a 



for T > 



(25) 



Here the 



background density in the vacuum is no(T, /x) = 0. For this model, we find that 
- + 1 — - = | and — = 1. Thus the critical exponent z = 2 and the correlation 

z uz 2 uz 1 

length exponent v = 1/2 with system dimension d = 1 can all be read off the universal 
scaling form (1251) . Meanwhile, the compressibility satisfies the scaling form 



k(T,h) - Ko(T,fi) w T«" 



•l- 



T 



(26) 



with Q(x) 



2e^ 



Li_i(— e x ) and Kq(T,ii) = 0. 



The equation of state ([2"T]) has been used to plot the universal low temperature 
quantum phase diagram in Figure 3. The vacuum state and the Tomonaga-Luttinger 
liquid persist below the two cross-over temperatures. They both vanish at the quantum 
critical point, as described by the scaling functions. Following Zhou and Ho |24j, 
the scaling functions for the density and compressibility can be used to map out the 
criticality of ultracold bosonic atoms in experiments. Figure 4 shows the detection of 
quantum criticality near \x c = using either the density or compressibility curves. We 
see clearly that such curves intersect at the critical potential /x c , indicative of scaling. 



4. Concluding remarks 

As highlighted in the Introduction, the motivation for the present work on one- 
dimensional bosons was the successful application of the polylog function method to 
the one-dimensional attractive Fermi gases [251 EE] of ultracold atoms up to order 1 /c 2 
[T4"t [T5] . Our key results are the equation of state (12T|) and the scaling forms (1231) and 
([23]) for the density and compressibility. The resulting low-temperature quantum phase 
diagram in Figure 3 is the simplest possible quantum criticality for quantum many-body 
systems. 

In general, the polylog function method is widely applicable to one- dimensional 
many-body systems with quadratic bare dispersion or linear bare dispersion e°(k) in 
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both attractive and repulsive regimes. The analytical polylog function method can 
thus play an important role in unifying the properties of attractive Fermi gases of 
ultracold atoms with higher symmetries. For example, the Yang- Yang method applied 
to the one- dimensional Fermi gas with attractive ^-function interaction and internal 
spin degrees of freedom leads to equations which may be reformulated according to the 
charge bound states and spin strings characterizing spin fluctuations [15] . For strong 
attraction, the spin fluctuations that couple to non spin-neutral charge bound states 
are exponentially small and can be asymptotically calculated [T4"| [15] . Thus the low 
energy physics is dominated by density fluctuations among the charge bound states. The 
full phase diagrams and thermodynamics of the one-dimensional attractive Fermi gases 
can be analytically calculated with relative ease using this polylog function method. 
For example, for spin- 1/2 attractive fermions, the universal Tomonaga-Luttinger liquid 
behaviour was identified from the pressure given in terms of polylog functions [14] . Here 
we further extend the results obtained using this approach. 

Following the procedure described in the previous section, we obtain the pressure 
p{T) = p b {T) + p u (T) from the TBA equations for spin- 1/2 attractive fermions [T3| 127]. 
The superscripts b and u denote bound pairs and excess fermions. For studying quantum 
criticality of ID attractive Fermi gas [IT], a high precision equation of state is highly 
desirable. To this end, the high order of corrections in T and l/\c\ are retained as 



p\T) 



T2 Li 3 -e Ab/T ) 





Li 3 (-e A "/ T ) 



+ 



(27) 



P U {T) = - 



-^TlLi3(-e^ T ) 
2nH 2 5 v > 



1 - 




+ o 4 ■ 



(28) 



Here we have defined the functions 



1 



tI 



2 c 4V2tt|c| 3 fh*\? - 



Ap u (T) 



4 T§ 



7T|C| 3 ( h 2 \§ 2 



\2mJ 

^U 5 (-e A °/ T ) 



(—)' 

\2m) 



Ti 



27r|c| 3 ( f? 



\2m) 



-Li,(-e^ T ) + / s 



(29) 



with 



A b =2fi + 



c 2 p\T) Ap u (T) 
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*._„ + £_*p +A . (30) 

The spin string contributions to the thermal fluctuations in the physically interesting 
regime, T < c 2 , T < H and (7! > 1, is given by f s = Te~"^e"~T^ I (7^), m 
which Jo is the standard Bessel function. 

The above results may be immediately applied to obtain accurate thermodynamic 
quantities such as the magnetization, specific heat and density profiles in a trapping 
potential. One can compare with and fit the experimental results obtained recently 
for trapped one-dimensional spin-1/2 fermions at Rice University j28|. Indeed, the 
experimental results confirm the expected phase diagram [221 EDI E3 ED E21 E31 EH E5] ■ 

On the other hand, for one-dimensional Fermi gases with repulsive interaction 
[35J EZ], we understand that antiferromagnetic effective spin-spin interaction directly 
couples to the charge degrees of freedom [38]. This triggers a spin-charge separated field 
theory of a Tomonaga-Luttinger liquid and an antiferromagnetic SU(N) Heisenberg 
spin chain. With the help of Wiener-Hopf techniques and polylog functions, we also 
can obtain a high precision equation of state for such one- dimensional repulsive Fermi 
gases [39]. The application of this approach to the study of thermodynamics of other 
one-dimensional degenerate Bose and Fermi gases, such as the spinor Bose gases and to 
mixtures of bosons and fermions, is also relatively straight forward. 

In general the equation of state provides essential insight into the thermodynamics 
of interacting many-body systems. Schemes have been proposed to directly measure 
the equation of state in experiments with ultracold atoms [UH HI]- Moreover, new 
schemes for mapping out the thermodynamics [42] and quantum criticality [21] of 
homogeneous systems by using the inhomogeneity of the trap can be directly applied to 
one-dimensional quantum many-body systems with a wide range of tunable interactions. 
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Interaction strength c 



Figure 1. Pressure of Lieb-Liniger bosons vs interaction strength c. The curves show 
the comparison between results obtained using the polylog function and numerical 
solution of the integral equation. 
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Temperature T Temperature T 



Figure 2. Finite temperature thermodynamics of Lieb-Liniger bosons. The curves 
show the comparison between the results obtained from numerics, the polylog function 
and Sommerfeld expansion for the entropy and the pressure. 
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Figure 3. Quantum phase diagram of Lieb-Liniger bosons. The plot shows the entropy 
in the T — U plane. The right red-dashed line denotes the crossover temperature 
separating the Tomonaga-Luttinger liquid (TLL) phase from the quantum critical 
regime (CR). The left red-dashed line separates the vacuum (V) from the quantum 
critical regime. For nonzero temperature, the vacuum produces particles with density 
n ~ je^^^ T with thermal wavelength A = ^raksT /2irh 2 which is much smaller than 
the inter-particle mean spacing. The vacuum can thus be taken as a semi-classical 
regime. 




Figure 4. Quantum criticality is mapped out from the density and compressibility 
at low temperatures. The left panel shows the compressibility vs chemical potential 
at temperatures T = O.OOleo, T — 0.003eo and T = 0.005eo- The inset shows the 
universal scaling function near the critical potential fi c = 0. The right panel shows 
intersecting density curves at the critical potential fi c = 0. 



